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ABSTRACT 

Short gamma-ray Bursts (SGRBs) are among the most luminous explosions in the universe, releasing in less 
than one second the energy emitted by our Galaxy over one year. Despite decades of observations, the nature 
of their "central engine" remains unknown. Considering a binary of magnetized neutron stars and solving the 
Einstein equations, we show that their merger results in a rapidly spinning black hole surrounded by a hot and 
highly magnetized torus. Lasting over 35 ms and much longer than previous simulations, our study reveals 
that magnetohydrodynamical instabilities amplify an initially turbulent magnetic field of ~ 10^^ G to produce 
an ordered poloidal field of ~ 10^^ G along the black-hole spin-axis, within a half-opening angle of ~ 30°, 
which may naturally launch a relativistic jet. The broad consistency of our ab-initio calculations with SGRB 
observations shows that the merger of magnetized neutron stars can provide the basic physical conditions for 
the central engine of SGRBs. 

Subject headings: Gamma-ray burst: general — black hole physics — stars: neutron — gravitational waves — 
magnetohydrodynamics (MHD) — methods: numerical 



1. INTRODUCTION 

The numerical investigation of the inspiral and merger of 
binary neutron stars (NSs) in full general relativity has made 
big strides in recent years. Crucial improvements in the 
formulation of the equations and numerical methods, along 
with increased computational resources, have extended the 
scope of early simulations. These developments have made 
it possible to compute the full evolution, from large binary- 
separations up to bl ack-hole (BH) formation, without and 
with mag netic field s (Shibata & Taniguchil2006UBaiotti et aLl 
2008 i And erson et a l. 2008 ; Liu et al. 2008 ; Giacomazzo et al.| 
2009, 2011 ), a nd with id ealized or realistic equations-of-state 
TeOS; Rezzolla et al.|2010i Kiuchi et al.,2010j . This tremen- 
dous progress is also providing information about the entire 
gravitational waveform, from the early inspiral up to t he ring - 
ing of the BH (see, e.g. [Duez| ( [2QT0l ); |Baiotti et al IpOTO] )). 
Advanced interferometric detectors starting from 201 4 are ex- 
pected to observe these sources at a rate of ~ 40 — 400 events 
per year ( Abadie et al. 2010). 

These simulations also probe whether the end-pr oduct of 
mergers can serve as the "central engine" o f SGRBs ( [Paczyn 



ski|1986tpchler et al.|1989]|Narayan et al.|1992j ). Thepreva- 
lent scenario invoked to explain SGRBs involves the coales- 
cence of a binary s ystem of compact objects, e.g. a BH and 
a NS or two NSs 
[20031 [NE||20Q7 



a NS or two NSs ([Ruffert & Jank a|[T999l [Rosswog et al. 

-Ruiz||2007|). After the 
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coalescence, the merged object is expected to collapse to a 
BH surrounded by an accretion torus. An essential ingredi- 
ent in this scenario is the formation of a central engine, which 
is required to launch a relativistic outflow with an energ y of 
^ 1q48 _ iq50 gj^g on a timescal e of ~ 0.1 — 1 s (N akar| 
[20071 |Lee & Ramir ez-Ruiz|[2QQ7]) With only one possiEIe 



exception (De Pasqu ale et al. 12 1(7), SGRB afterglows do not 
clearly show a jet- associated light-curve steepening ( [Nakar 
[2007 ), thus suggesting typical jet opening half-angles of at 
least several degrees. 

The qualitative scenario described above is generally sup- 
ported by the association of SGRBs with old stellar popula- 
tions, distin ct from the young massive star ass ociations for 
long GRBs ( |Fox et al.|2005l |Prochaska et al.^2006). It is also 
supported to a good extent by fully general-relativistic sim- 
ulations, which show that the formation of a torus of mass 
Mtor < 0.4 around a BH with spin J/M^ - 0.7 - 0.8, 
is inevitable ( [Rezzolla et al.|2010| ). However, the simulations 
have so far failed to show that a jet can be produced. We 
here provide the first evidence that the merger of a binary of 
modestly magnetized NSs naturally forms many of the con- 
ditions needed to produce a jet of ultrastrong magnetic field, 
with properties that are broadly consistent with SGRB obser- 
vations. 

2. NUMERICAL SIMULATIONS 

For the simulations we use the Cactus/CarpetAVhisky codes 
(^Schnetter et al. [[20 02 : Thornburg 2004 ; Giacomazzo e t al.[ 
,2007; PoUney et al. |200 7 ) and we consider a configuration 
that could represent the properties of a NS -binary a few or- 
bits before its coalescence, within a fully general-relativistic 
and an ideal-ma gnetohydrodynamic (MHD) framework | Gi-[ 
acomazzo et al. 2007] [2011|. More specifically, we simu- 
late two equal-mass NSs, each with a gravitational mass of 
1.5 {i.e. sufficiently large to produce a BH soon after 
the merger), an equatorial radius of 13.6 km, and on a cir- 
cular orbit with initial separation of 45 km between the 
centres (all lengthscales are coordinate scales; (Taniguc hi &[ 
Gourgoulhon 2002) ). Confined in each star is a poloidal mag- 
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Figure 1. Snapshots at representative times of the evolution of the binary and of the formation of a large-scale ordered magnetic field. Shown with a color-code 
map is the density, over which the magnetic-field lines are superposed. The panels in the upper row refer to the binary during the merger (t = 7A ms) and before 
the collapse to BH (t = 13.8 ms), while those in the lower row to the evolution after the formation of the BH (t = 15.26 ms, t = 26.5 ms). Green lines sample 
the magnetic field in the torus and on the equatorial plane, while white lines show the magnetic field outside the torus and near the BH spin axis. The inner/outer 
part of the torus has a size of ~ 90/170 km, while the horizon has a diameter of 2:^ 9 km. 



netic field with a maximum strength of 10 ^^ G (indicated as 
Ml . 62-B12 in G iacomazzo et al.||2QTT] ). At this separa- 
tion, the binary loses energy and angular momentum via emis- 
sion of gravitational waves (GWs), thus rapidly proceeding on 
tighter orbits as it evolves. After about 8 ms (~ 3 orbits) the 
two NSs merge forming a hypermassive NS (HMNS), namely, 
a rapidly and differentially -rotating NS, whose mass, 3.0 M©, 
is above the maximum mass, 2.1 M©, allowed with uniform 
rotation by our ideal-gas EOS^ with an adiabatic index of 2. 
Being metastable, a HMNS can exist as long as it is able to 
resist against collapse via a suitable redistribution of angu- 
lar mome ntum (e.g. deforming into a "bar" shape, Shibata &] 
[Taniguchi 2006; B aiotti et al.||2008| ), or through the pressure 
support coming from the large temperature-increase produced 
by the merger. However, because the HMNS is also losing an- 
gular momentum through GWs, its lifetime is limited to a few 
ms, after which it collapses to a BH with mass M = 2.91 Mq 
and spin J/M^ = 0.81, surr ounded by a hot and dens e torus 
with mass M^or = 0.063 ( [Giacomazzo et al.pOll] ). 



3. DYNAMICS OF MATTER AND MAGNETIC FIELDS 

These stages of the evolution can be seen in Figure [T] 
which shows snapshots of the density color-coded between 
10^ and 10^^ gr/cm^, and of the magnetic field lines (green 

^ The use of a simplified EOS does not influence particularly our results 
besides determining the precise time when the HMNS collapses to a BH. 



on the equatorial plane and white outside the torus). Soon 
after the BH formation the torus reaches a quasi- stationary 
regime, during which the density has maximum values of 
~ lO^^g/cm^, while the accretion rate settles to M c:^ 
0.2Mq/s. Using the measured values of the torus mass 
and of the accretion rate, and assuming the latter will not 
change significantly, such a regime could last for taccr = 
Mtor/M 0.3 s, after which the torus is fully accreted; fur- 
thermore, if the two NSs have unequal masses, tidal tails are 
produced w hich provide additional late-time accretion ( .Rez-| 



zolla et aiypO lO). This accreti on timescale is close to the typi- 
cal ob serves SGRB durations ( [Kouveliotou et al.|1993t|Nak^ 



2007| ). It is also long enough for the neutrinos produced in 
the torus to escape and annihilate in its neighborhood; es- 
timates of the associated energv deposition rate range from 
^ 10^^ erg/s (D essart et al.|2QQ9| ) to ^ 10^° erg/s ( Setiawan] 
et al. 12004,), thus leading to a total energy deposition between 



aTewTO^erg and a few 10^^ erg over a fraction of a second. 
This energy would be sufficient to launch a relativistic fire- 
ball, but because we do not yet account for radiative losses, 
the large reservoir of thermal energy in the torus cannot be 
extracted in our simulations. 

The GW signal of the whole process is shown in the left 
panel of Figure |2] while the bottom part exhibits the evolu- 
tion of the MHD luminosity, Lmhd, as computed from the 
integrated Poynting flux (solid line) and of the correspond- 
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Figure 2. Left panel: GW signal shown through the ^ = 2, m = 2 mode of the + polarization, (/z+)22, (top part) and the MHD luminosity, Lmhd. (bottom 
part) as computed from the integrated Poynting flux and shown with a solid line. The corresponding energy, £^mhd. is shown with a dashed line. The dotted and 
dashed vertical lines show the times of merger (as deduced from the first peak in the evolution of amplitude of the GW amplitude) and BH formation, respectively. 
Right panel: Evolution of the maximum of the magnetic field in its poloidal (red solid line) and toroidal (blue dashed line) components. The bottom panel shows 
the maximum local fluid energy indicating that an unbound outflow (i.e. E\oc > 1) develops and is sustained after BH formation. 



ing energy, Emyld (dashed Hne). Clearly, the MHD emission 
starts only at the time of merger and increases exponentially 
after BH formation, when the GW signal essentially shuts 
off. Assuming that the quasi-stationary MHD luminosity is 
4 X 10^^ erg/s, the total MHD energy released during the 
lifetime of the torus is 1.2 x 10^^ erg, which, if spread over 
an opening half-angle of ~ 30° (see discussion below), sug- 
gests a lower limit to the isotropic equivalent energy in the 
outflow of :^ 9 X 10^^ erg. While this is at the low end of 
the observed distribution of gamma-ray energies for SGRBs, 
larger MHD luminosities are expected either through the ad- 
ditional growth of the magnetic field via the ongoing winding 
of the field lines in the disk (the simulation covers only one 
tenth of taccr), or when magnetic reconnection (which cannot 
take place within our ideal-MHD approach), is also accounted 
for (which may also increase the gamma-ray efficiency; see, 
e.g. McKinney & Uzdensky 2010). 

The last two panels of Figure [T] offer views of the accreting 
torus after the BH formation. Although the matter dynamics 
is quasi- stationary, the last two panels clearly show that the 
magnetic field is not and instead evolves significantly. It is 
only when the system is followed well after the formation of 
a BH, that MHD instabilities develop and generate the cen- 
tral, low-density, poloidal-field funnel. T his regime, which 
was not accessible to previous simulations ([Price & Rosswog 



[20061 [Anderson et al.|2008t|Liu et al.||2008|), is ess ential tor 
the jet formation ( Aloy et al.[[2005 Komissarov et al. 2009). 
Because the strongly magnetized matter in the torus is highly 
conductive, it shears the magnetic-field lines via differential 
rotation. A measurement of the angular velocity in the torus 
indicates that it is essentially Keplerian and t hus unstable to 
the m agneto-rotational instability (MRI; Bal bus & Hawley 
1998|), which develops 5 ms after BH formation and ampli 



ties exponentially both the poloidal and the toroidal magnetic 
fields; the ^-folding time of the instability is 2.5 ms and in 
good agre ement with the one ex pected in the outer parts of the 
torus (Bal bus & Hawley[|1998[ ). Because of this exponential 
growth, the tinal value or the magnetic field is largely insen- 
sitive to the initial strength and thus a robust feature of the 



dynamics. 

A quantitative view of the magnetic-field growth is shown 
in the right panel of Figure [2j which shows the evolution 
of the maximum values in the poloidal and toroidal com- 
ponents. Note that the latter is negligibly small before the 
merger, reaches equipartition with the poloidal field as a result 
of a Kelvin-Helmholtz instabil ity triggered by the shearing of 
the stellar surfaces at merger ( Price & Rosswog|2006[[Giaco^ 



mazzo et al.|2009| , and finally grows to 10^^ G by the end 



of the simulation. At later times {t > 22 ms), when the insta 
bility is suppressed, the further growth of the field is due to 
the shearing of the field lines and it increases only as a power 
law with exponent 3.5 (4.5) for the poloidal (toroidal) com- 
ponent. Although the magnetic-field growth essentiall y stalls 
after t 35 ms, fur ther slower growths are possible ( [Qber- 



'gaulin ger et al.|200"9] ), yielding correspondingly larger Poynt- 
ing fluxes. Indeed, when the ratio between the magnetic flux 
across the horizon and the mass accretion rate becomes suf- 
ficiently large , a Blandford-Zn ajek mechanism (Blandford &| 
[Znajek[[1977| ) may be ignited ( [Komissarov & Barkov[[2009'] ); 
such conditions are not met over the timescale of our simula- 
tions, but could develop over longer timescales. Also shown 
in the right panel of Figure |2] is the maximum local fluid en- 
ergy, highlighting that an unbound outflow (i.e. Ei^c > 1) de- 
velops after BH formation along the outer walls of the torus 
and persists for the whole duration of the simulation. 

Finally, Figure [3] provides a summary of the magnetic-field 
dynamics. It shows the magnetic field in the HMNS formed 
after the merger and its structure and dynamics after the col- 
lapse to BH. In particular, in the last three panels it shows 
the magnetic-field structure inside the torus and on the equa- 
torial plane (green), and outside the torus and near the axis 
(white). It is apparent that the highly turbulent magnetic field 
in the HMNS (t = 13.8 ms) changes systematically as the 
BH is produced (t = 15.26 ms), leading to the formation 
of a toroidal magnetic field in the torus (t = 21.2 ms). As 
the MRI sets in, the magnetic field is not only amplified, but 
also organizes itself into a dual structure, which is mostly 
toroidal in the accretion torus with ^tor — 2 x lO^^G, but 




Figure 3. Magnetic-field structure in the HMNS (first panel) and after the collapse to BH (last three panels). Green refers to magnetic-field lines inside the 
torus and on the equatorial plane, while white refers to magnetic-field lines outside the torus and near the axis. The highly turbulent, predominantly poloidal 
magnetic-field structure in the HMNS (t = 13.8 ms) changes systematically as the BH is produced (t = 15.26 ms), leading to the formation of a predominantly 
toroidal magnetic field in the torus (t = 21.2 ms). All panels have the same linear scale, with the horizon diameter being of :^ 9 km. 



predominantly poloidal and jet-like along the BH spin axis, 
with 5poi 8 X 10^"^ G (t = 26.5 ms). Note that the gen- 
eration of an ordered large-scale field is far from trivial and 
a nonlinear dynamo may explain why the MRI brings a mag- 
netic field self-organization, as it has been also suggested in 
case of MRI-mediated growth of the magnetic field i n the con- 
ditions met in the collapse of massive ste llar cores ( |Lesur &| 
Qgilvie|[2QQ8l |Qbergaulinger et al.||2QQ9| ). However, the jet- 
like structure produced in the simulation is not yet the highly 
collimated ultrarelativistic outflow expected in SGRBs (see 
also below). 

The hollow jet-like magnetic structure has an opening half- 
angle of ^ 30°, which sets an upper limit for the opening half- 
angle of any potential outflow, either produced by neutrino 



energy deposition ( Alo y et ar]|2Q05|) or by electromagnetic 
(EM) processes ( Komis sarov et al . |2009| ) . In our simulations 



most of the outflow develops along the edges of the jet-like 
structure, via a turbulent layer of EM driven matter, which 
shields the central funnel from excessive baryonic pollution. 
We envision that such a layer is crucial to set the opening 
angle of any ultrarelativistic jet, to shape both the radial and 
transverse structure of the jet, as well as to determine its sta- 
bility properties. The Lorentz factors of the outflow measured 
in our simulations are not very high (F < 4), but can poten- 
tially be amplified by several orders of magnitude in the inner 
baryon-poor re gions through special-relativistic effects ( Aloy 
& Rezzolla|200 6) or the variability of the flow (GranotetaLj 
201 We expect that such accelerations will be produced as 



a more realistic and general-relativistic treatment of the radia- 
tive losses will become computationally affordable. 

4. COMPARISON WITH OBSERVATIONS 

Below we briefly discuss how our results broadly match the 
properties of the central engine as deduced from the observa- 
tions. 



4.1. Duration 

The observed duration of the prompt gamma-ray emission 
GRBs is energy dependent and is usually determined through 
Tx, the time over which x% of the total counts are observed, 
between the (100 — x)/2 and (100 + x)/2 percentiles. The 
most common intervals used are Tgo (or T50), initially de- 
fined (KouveHotou et al. 1993) between 20 keV and 2 MeV. 
The G RB duration distribution is bimodal (Kouveliotou etaL] 
119931), where the durations of SGRBs (approximately 25% of 
GRBs) are well-fit by a fairly wide log-normal distribution 
centered around T90 ~ 0.8 s with a FWHM of 1.4 dex ( |Nakar| 
2007) ). The typical redshifts of the SGRBs observed with 
Swift are in the range 2: ~ 0.3— 1, suggesting a central value of 
the intrinsic duration distribution of^ (l+z)~^0.8s~0.5s, 
and a comparably wide distribution around this value. This is 
in close agreement with our accretion time of ~ 0.3 s. 



4.2. Energy 
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The isotropic equivalent energy output in the prompt 
gamma-ray emission of SGRBs, E^^i^o^ spans a wide range, 
from (2.7 =b 1) x 10^^ erg (in the observed energy range 
15 - 350 keV) for ORB 050509B at a redshift of z = 
0.225 ( [Bloom et al.l|2006| ), up to (1.08 ± 0.06) x 10^^ erg 
(in the observed energy ra nge 10 keV - 30 G eV) for 
GRB 090510 at ^ = 0.903 ( |Ackermann eraL||2QTQl ). How- 
ever, the most typical va lues are in the range ^^,iso ^ 10^^ — 
10^1 erg (Nakarl 2007). In our model the highly relativis- 
tic outflow may be powered either by neutrino-anti neutrino 
annihilation, or by the Blandford-Znajek mechanism. For 
the former one might expect a t otal energy relea se between 
a few 10"^^ erg and - 10"^^ erg ( [Oechslin et aLl|2006[ [Birkl 



et al.|200 7 ), into a bipolar relativistic jet of opening half- angle 
^jet ^ 8 — 30°, corresponding to a fraction /5~0.01 — 0.13 
of the total solid angle, and isotropic equivalent energies, 
Evy;iso^ between a few 10^^ erg and ~ 10^^ erg. For the 
latter mechanism, instead, and if the magnetization near the 
event horizon becomes sufficiently high, we could expect a jet 
power for our values fo r the BH mass and spin that is of ( [Lee 
let al.|2000, Perez-Ramirez et al.|2010| ) 



Lbz ^ 3.0 X 10 



50 



/rel 
0.1 



B 



2 X 1015 G 



erg/s, (1) 



where /rei is the fraction of the total Blandford-Znajek power 
that is channeled into the resulting relativistic jet (and /rei ~ 
0.1 might be expected for ejecta with asymptotic Lorentz fac- 
tors above 100). This relativistic outflow is launched over a 
timescale of ~ 0.2 s and corresponds to 
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51 
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B 



2 X 1015 G 



erg . 
(2) 

Comparing the X-ray afterglow luminosity (after 10 or 11 
hours) and £^^,iso suggests that the efficiency o f the prompt 
gamma-ray emissio n in SG RBs is typically high dBloo m et"aL] 
[2006; Gehrels et al. |2009[ ), and similar to that of long GRBs 
jGranot et al., 2006), with £;^,iso ^ (0.1 - 0.9)£;iso, radiat- 
ing between ^ 10% and ~ 90% of the initial energy of the 
ultrarelativistic outflow. Therefore, our model is able to ac- 
commodate the observed £^^,iso values. 

4.3. Lorentz factor 

The Fermi Gamma-R ay Space Telescope has detected GeV 
emission from SGRBs ( jAbdo et al.||2009a| ), suggesting typi- 
cal lower limits of Fmin ^ 10"^ — 1 0"^ . In particular, Fmin ^ 
1200 was obtained for GRB 09051 ([Ackermann et al |2QTo1 ). 
However, a more realistic model ( [Granot et al.|2008|) results 
in Fmin values lower by a factor of ^ 3 dAcker mann et al.j 
2011] ). Therefore, the central engine should be capable of 
producing outflow Lorentz factors of at least a few hundred. 
The fact that our simulation produces a strongly magnetized 
mildly relativistic outflow at angles near ~ 30° from the BH 
spin axis would help shield the inner region near the spin axis 
from excessive baryon loading, and thus assist in achieving 
high asymptotic Lorentz factors at large distance from the 
source, after the outflow in this region is triggered by neu- 
trinos and/or the Blandford-Znajek mechanism. 

4.4. Jet angular structure 

This is poorly constrained by observations (even more so 
than for long GRBs). The only compelling case for a jet 



break in the afterglow light curve is for GRB 090510 ( |De| 
Pasquale et al.[[20TQ| ), which occurred very early on (after 
~ 1400 s), and would thus imply an extremely narrow jet 
(^jet 0.2 — 0.4°) and modest true energy output in gamma 
rays (~ 10^^ erg). If this is indeed a jet break, it might cor- 
respond to a line of sight near a very narrow and bright core 
of a jet, which also has significantly wider wings. Observers 
with lines of sight along these win gs would then see a much 
diminer an d more typical SGRB ( [Rossi et"aL| |2002[ [Peng 
et al.| [2005'i Tlacusin et al. 2009); without such wings, how- 
ever, tRe observations would suggest a very large intrinsic and 
beaming-corrected event rate per unit volume. In mos t cases 
there are only lower limits on a possible jet break time (jNakar 



2007), resulting in typical limits of > 10 ^ or 



jet 



1 his is consistent with our expectation of 6>jet ^ 8 — 30° 
for the ultrarelativistic ejecta capable of producing a SGRB 
(which would also imply a reasonable SGRB intrinsic event 
rate per unit volume). 

5. CONCLUSIONS 

The calculations reported here clearly demonstrate that a 
binary merger of two NSs inevitably leads to the formation 
of a relativistic jet-like and ultrastrong magnetic field, which 
could serve as a central engine for SGRBs. Because the 
magnetic-field growth is exponential, the picture emerging 
from our simulations is rather general and applies equally 
even to mildly magnetized NSs. Overall, this work removes a 
significant earlier uncertainty as to whether such binary merg- 
ers can indeed produce the central engines of SGRBs. While 
the EM energy release is already broadly compatible with the 
observations, our simulations lack a proper treatment of the 
energy losses via photons and neutrinos, which can provide 
a fundamental contribution to the ener gy input necessary to 
launch the fi reball and cool the torus ( Setiawan et al.| [2004t 
Dessart et al . 2009 ). This additional energy input, whose self- 
consistent inclusion in general relativity remains extremely 
challenging, may help to launch an ultrarelativistic outflow 
very early after the BH forms and complete the picture of the 
central engine of a SGRB. 
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